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Surface modes have been studied in different branches of physics including guided 
wave optics, where surface waves were predicted to exist at interfaces separating 
periodic and homogeneous dielectric mediaP Recently, it was suggested theoreticallj'^ 
and demonstrated experimentalljJ^ that nonlinearity-mediated trapping of light near 
the edge of a truncated waveguide array with self-focusing nonlinear response can 
lead to the formation of nonlinear localized surface states which can be understood 
as discrete optical soliton^ localized near and trapped by the surfac^ for powers 
exceeding a certain threshold value. One of the important generalizations of these 
ideas is the concept of multi-gap surface solitons, i.e. mutually trapped surface states 
with components associated with different spectral gaps.l^ 

In this Letter, we study another important generalization of the concept of non- 
linear surface modes. We analyze linear and nonlinear optical guided modes localized 
at an interface separating two different semi-infinite periodic photonic lattices. In the 
framework of an effective discrete model we demonstrate that the analysis of linear 
interface states in such composite arrays provides an important tool for analyzing the 
interface solitons and their basic properties. We then find numerically the families of 
stable interface lattice solitons including a novel class of hybrid staggered/unstaggered 
lattice solitons with tails localized in spectral gaps of different types. 

We consider an interface separating two different semi-infinite arrays of optical 
waveguides (as shown in the top of Fig. 1) described by the system of coup led- mode 
equation^ for the normalized mode amplitudes En- 

dE 

i—^ + e„E„ + + K-i) + l\En?E^ = 0, (1) 

dz 
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where the propagation coordinate z is normahzed to the intersite couphng V ^ are 
defined in terms of the actual electric field Sn as En — {2V Xqttiq / T:nQn2Y^'^ Sn, Aq is the 
free-space wavelength, 770 is the free-space impedance, no and n2 are the mean values 
of the linear and nonlinear refractive indices of each waveguide, and 7 defines the 
nonlinear response strength. The waveguide interface is introduced by the condition: 
eo at n = 0, and e„ = or e„ = for negative or positive n, respectively. 

First, we look for linear surface modes in the form E^ — A ^1^' exp{if3z) localized 
near the interface waveguide with n — 0, and obtain the condition = e^o/eBo 

and the dispersion relation (3 — + |(e^o + cbo)(1 ~ y^l -|- 4/eAoeBo), where e^o = 
ea — Co and e^o = es — Co- Figure 1 summarizes our results for the existence of such 
localized states on the parameter plane (e^o; ^bo), ^ well as displays examples of 
localized modes corresponding to different existence regions. We note the existence of 
two sectors where no locahzed states exist (shaded regions) . One of them is bounded 
by the curves sbo = 4/|eAo| (for -00 < e^o < 0) and sbo = eAo/(l - cao) (for 
—00 < eao < !)• The other one is obtained by a reflection across the origin. Inside 
these regions either or both exceed one. 

One of the important observations that follows from our analysis is the existence 
of hybrid staggered/unstaggered interface modes for the opposite sign of the propaga- 
tion constant mismatches of two lattices. These modes have the tails localized in the 
bandgaps of different types, i.e. above (for one array) and below (for the other array) 
of the spectral band. 

In Figs. 2(a,b) we show the propagation constant /? of the localized modes as a 
function of the interface parameter cq, for the characteristic cases (a) = 0.6, = 
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Fig. 1. (Color online) Phase diagram of different types of localized interface 
modes. No localized modes exist inside the shaded areas. The insets show 
examples of localized modes corresponding to different values of bao = — co 
and eBo = — eo- Top: schematic structure of the waveguide interface. 

—0.6, and (b) = 3, = —3. We note that the mode always lies outside the linear 
spectral bands, but its structure depends strongly on the overlap of the bands, so that 
the hybrid modes appear for a relatively large band mismatch, as shown in Fig. 2(b) 
(middle curve). 

The analysis of linear localized interface modes in such an array provides an im- 
portant information about the existence of nonlinear interface modes — lattice surface 
solitons. Next, we consider two semi-infinite nonlinear waveguide arrays character- 
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Fig. 2. (Color online) Families of localized modes shown as the dependencies 
(3 vs. eg, for the cases: (a) = 0.6, = —0.6, and (b) = 3,6^ = —3. No 
localized modes exist inside the shaded regions. The dashed lines mark the 
spectral bands of both arrays, which in the case (b) do not overlap with each 
other. 



ized by propagation constants and that are joined by an interface waveguide 
with the propagation constant eo. We focus on the interface modes defined by having 
their centers at either the first of the A waveguides or the first of the B waveguides. 
We find different classes of nonlinear localized modes and study their linear stability 
numerically via a multidimensional Newton-Raphson method. 

First, we consider the case eo = = —^b = 0.6 and 7 = +1. In the linear limit, 
this corresponds to a negative vertical line in Fig. 1, where no localized modes exist. 
The presence of nonlinearity shifts the propagation constant of the mode to the left, 
until it gives rise to an unstaggered localized mode. Therefore, we predict the lowest 
nonlinear interface mode to be unstaggered. This is indeed confirmed by our numerical 
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Fig. 3. (Color online) Power vs. propagation constant for the interface unstag- 
gered solitons centered on the first A waveguide (for = eo = 0.6, = —0.6). 
The solid (dashed) curve refers to the stable (unstable) branches. Inserts show 
two examples of the interface modes. 



computations, and the family of the lowest-order interface nonlinear modes is shown 
in Fig. 3, where the upper/lower branch corresponds to unstable/stable modes. Next, 
we consider the case = 3, eo = and = —3. Results are summarized in Fig. 4 
which shows the dependence of the power vs. propagation constant for several low- 
power modes. We note that the lowest mode extends all the way to zero power, and 
therefore it corresponds in that limit to the linear mode induced by three concurrent 
dissimilar propagation constants. More importantly, the mode amplitudes show now 
a hybrid character, being unstaggered in one side of the interface and staggered on 
the other. 

In addition, we find many other types of nonlinear interface modes including 
twisted and fiat-top modes, as well as the modes with different location of their 
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Fig. 4. (Color online) Power vs. propagation constant for the hybrid inter- 
face staggered/unstaggered lattice solitons for eA = — 3, eo = 0,6^ = 3. Solid 
(dashed) curves refer to the stable (unstable) branches. Inserts show three 
examples of the hybrid interface solitons. 

centers relative to the interface, as discussed earlier for a semi-infinite waveguide 
array.l^ As a special limit of those modes, we find also the kink surface modes which 
are extended in one direction while being localized in the other. 

Next, we study the evolution of all types of nonlinear interface modes checking 
in this way our stability results and analyzing the evolution scenario for the unstable 
modes. In all cases examined, we observe that the unstable modes decay into the 
unstaggered fundamental 

mode by emitting radiation. In particular, for the unstaggered unstable modes, 
this decay is rapid and it proceeds vertically. For other cases, the unstable localized 
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Fig. 5. (Color online) Top: trapping efficiency for the generation of interface 
solitons. The red (blue) curve denotes the case when the initial input is cen- 
tered on the ffist A (B) waveguide. Bottom: evolution of initial states marked 
'a' and 'b' (for = 0.6, = -0.6). 



modes decay into higher-power generalizations of the fundamental mode. Also, we 
observe that unstable unstaggered modes decay much faster than unstable twisted, 
flat-top, or dark-like nonlinear modes. 

Finally, we study how an interface lattice soliton can be generated in experiment 
by exciting a single waveguide which is either the very first of the A waveguides, or the 
very first of the B waveguides. We define the trapping efficiency as the power fraction 
Pout /-Pin remaining in the 10 central waveguides. Results for the trapping efficiency 
are shown in Fig. 5(top), where we have used 201 waveguides, a total evolution length 
of 20. The bottom portion of Fig. 5 shows the evolution of the initial states marked 
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by the points 'a' and 'b' in Fig. 5(top). 

We would like to emphasize that the results obtained here can be easily gen- 
eralized to the case of surface gap solitons which were predicted theoreticalljP' and 
observed experiment all j0 in periodic photonic lattices with defocusing nonlinearity, 
where surface solitons appear in the gaps of the photonic bandgap spectra or their 
overlaps. 

In conclusion, we have analyzed different types of linear and nonlinear optical 
guided modes localized at the interface separating two different semi-infinite periodic 
photonic lattices. In the framework of an effective discrete nonlinear model, we have 
demonstrated the existence of stable interface lattice solitons including the hybrid 
staggered/unstaggered discrete solitons with the tails that belong to different spectral 
gaps. We believe our results will encourage the first experimental observations of this 
novel type of surface optical solitons in photonic lattices. 

This work has been supported by Conicyt and Fondecyt grants 1050193 and 
7050173 in Chile, and by the Australian Research Council in Australia. 
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List of Figure Captions 

Figure 1: (Color online) Phase diagram of different types of localized interface modes. 
No localized modes exist inside the shaded areas. The insets show examples of local- 
ized modes corresponding to different values of e^o = ^a — and e^o = — ^o- Top: 
schematic structure of the waveguide interface. 

Figure 2: (Color online) Famihes of localized modes shown as the dependencies P 
vs. eo, for the cases: (a) €a — 0.6, — —0.6, and (b) — 3,eB — —3. No localized 
modes exist inside the shaded regions. The dashed hues mark the spectral bands of 
both arrays, which in the case (b) do not overlap with each other. 

Figure 3: (Color onhne) Power vs. propagation constant for the interface unstaggered 
solitons centered on the first A waveguide (for ca — eo — 0.6, = —0.6). The solid 

(dashed) curve refers to the stable (unstable) branches. Inserts show two examples of 
the interface modes. 

Figure 4: (Color online) Power vs. propagation constant for the hybrid interface 
staggered/unstaggered lattice solitons for €a — — 3, eo = 0,6^ = 3. Solid (dashed) 
curves refer to the stable (unstable) branches. Inserts show three examples of the 
hybrid interface solitons. 

Figure 5: (Color online) Top: trapping efficiency for the generation of interface soli- 
tons. The red (blue) curve denotes the case when the initial input is centered on 
the first A (B) waveguide. Bottom: evolution of initial states marked 'a' and 'b' (for 
ca — 0.6, cb — —0.6). 
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